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Dynamical decoupling protocols are one of the most used tools for efficient quantum error correc¬ 
tions and for reservoir engineering. In this paper we study the effect of dynamical decoupling pulses 
on the preservation of both quantum and classical correlations, and their influence on the intriguing 
phenomenon of time-invariant discord. We study two qubits experiencing local dephasing with an 
Ohmic class spectrum and subject to dynamical decoupling protocols. We investigate the connection 
between the dynamics of both classical and quantum correlations and the behaviour of information 
flow between system and environment. Finally, we establish a set of necessary conditions for which 
time-invariant discord can be created using dynamically decoupling techniques. 
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I. INTRODUCTION 

Genuine quantum traits possessed by quantum me¬ 
chanical systems, such as coherence and quantum cor¬ 
relations, are fundamental for all possible applications 
of quantum information science. Realistic quantum sys¬ 
tems are however subjected to the interaction with their 
surroundings, which consistently results in the loss of 
such critical quantum features [T] . Since quantum devices 
need to be robust against the destructive effects of their 
environment to be used in everyday tasks, protecting co¬ 
herence and quantum correlations as long as possible is 
both of great importance and a major challenge. 

Several methods have been introduced to suppress the 
noise induced on the system by the environment in or¬ 
der to preserve quantum correlations. One such strat¬ 
egy, known in the literature as dynamical decoupling 
(DD), relies on the application of a sequence of exter¬ 
nal pulses to the principle system so that the harmful 
effects of the environment are countered H? -ig. DD 
techniques for open quantum systems are among the 
most successful methods to subdue decoherence. Specifi¬ 
cally, “bang bang” periodic dynamical decoupling (FDD) 
schemes can be exploited to prolong coherence times and 
restore monotonically decaying correlations in quantum 
systems which are undergoing decoherence HIT]. 

Another approach to counteract the environmental 
noise is based on the exploitation of memory effects 
characterising non-Markovian dynamics While in 

Markovian open systems information is monotonically 
lost to the environment, non-Markovian dynamics ex¬ 
hibits backflow of information from the environment to 
the system. This in turn results in the emergence of mem¬ 
ory effects since, as a consequence of the backflow, future 
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states of the system might depend on its past states. 
Thus, similarly to DD techniques, non-Markovianity can 
also be used to restore monotonically decaying correla¬ 
tions m- An ongoing debate concerns itself with char¬ 
acterizing and quantifying non-Markovianity |11H19) . At 
present, no universal measure exists and indeed in some 
cases, connections between measures can not be estab¬ 
lished EOIIII]. However, our treatment here is not quan¬ 
titatively dependent on a specihc measure, allowing our 
results to be presented independently of a specific quan- 
tiher of non-Markovianity. 

Entanglement has been considered as the main re¬ 
source of quantum protocols for a long time [22]. How¬ 
ever, it has been later on demonstrated via numerous 
examples that quantum discord |23| is also essential for, 
e.g., distribution of entanglement |24j . quantum lock¬ 
ing [25] , entanglement irreversibility [26] and many other 
tasks m- Quantum discord has been shown to be frozen 
at a finite value for a certain time interval |28] or even 
forever [33], depending on the properties of the initial 
state and the specifics of the interaction with the envi¬ 
ronment. In particular, for two qubits in Bell-diagonal 
states and independently interacting with dephasing en¬ 
vironments, while frozen discord manifests itself for both 
Markovian [28] and non-Markovian [30] environments, 
the emergence of time-invariant discord is more closely 
related to the presence of memory [33] . 

The use of either DD techniques or non-Markovian 
effects prolongs the preservation of both entanglement 
and discord in presence of environmental noise [3TVI35] . 
However, studies of single qubit dynamics have recently 
shown that the simultaneous use of non-Markovian reser¬ 
voir engineering and DD protocols is counterproductive 
for coherence preservation [7]. Motivated by this rather 
unexpected result, here we study the possibility of cre¬ 
ating time-invariant discord through FDD techniques in 
relation to the non-Markovian behavior in the dynam¬ 
ics. A clear merit in creating time-invariant discord is 
in the execution of quantum protocols based on discord, 
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which become unaffected by noise under the specific con¬ 
ditions. In this work, we consider two qubits interact¬ 
ing with a dephasing environment having an ohmic type 
spectral density, focussing on the cases of both one-sided 
and two-sided local noise. Our investigation allows us 
to establish a set of necessary conditions for which time- 
invariant discord can be observed. 

The paper is structured as follows. Section II in¬ 
troduces the system, namely a pure dephasing model, 
along with its exact solution in presence of DD. In Sec¬ 
tion III, we state the main results of our study, i.e., 
we present a detailed investigation of the phenomenon 
of time-invariant discord and determine a set of neces¬ 
sary conditions for which time-invariant discord mani¬ 
fests when the system is periodically pulsed. In Section 
IV we present conclusions. 

II. THE MODEL 

We consider a microscopic model describing the de¬ 
phasing interaction of two identical qubits with their re¬ 
spective local independent environments represented by 
a bath of bosonic oscillators. The total Hamiltonian of 
each of the two independent qubit-environment parts can 
be written as follows 

H = ^ uika\ak + ^ cr^ighak + ( 1 ) 

k k 

in units of h. Here, is the z-component of the spin op¬ 
erator for the two-level systems, with wq the qubit transi¬ 
tion frequency and uik the frequency of the fcth reservoir 
mode. The corresponding annihilation (creation) opera¬ 
tors are given by afc(a^), and gk is the coupling constant 
associated with the qubit-fcth mode interaction. This 
model is exact and admits an analytical solution [36II38j . 
In the interaction picture, the master equation for the 
density matrix of the qubit is given by 

P = loit)[(J^pcr^- p]/2, (2) 

The resulting decay of the off-diagonal density matrix 
elements are as follows: poi{t) = Pio{t) = poi{0)e~^°^*^ 
where 70 (t) = dTo{t)/dt and 

^oit)= [ ^^^[I - cos(w<)]dw, (3) 

Jo ^ 

where /(w) = J2j ~ the spectral density 

function characterising the interaction of the qubit with 
the oscillator bath, which is assumed to be initially at 
zero temperature. In our work, we consider the widely 
studied class of spectral densities of the form 

(4) 

UJc 

with s being the Ohmicity parameter and uJc a cutoff 
frequency. Ohmic spectra correspond to the case s = 1, 
sub-Ohmic to s < I and super-Ohmic to s > I. 


Next, we introduce the analytical solution of the above 
dephasing model in the presence of an arbitrary se¬ 
quence of bang-bang pulses modelled as instantaneous tt- 
rotations about the x-axis. In this situation, the reduced 
dynamics of the qubit can still be exactly described by 
exchanging ro(<) with a modified decoherence function 
r(t). Suppose that we consider an arbitrary storage time, 
t, during which a total number of N pulses are applied at 
instants {ti, ...t/}, with 0 < ti < t 2 < ■■■ < tf < t. 
Then, the decoherence function r(t) can be expressed as 

m 

( Toit) t < 

r(t) = < r„(t) tn < t < tn+ 1,0 < n < N , ( 5 ) 

[ rAr(t) tf <t 

where the exact representation of the controlled decoher¬ 
ence function r„(t) for 1 < n < iV, can be written in the 
following form: 

roo Tl^.\ 

^n{t) = J -^^\yn{ujt)\'^duj, n>0, (6) 

where |yo(wt)P = |1 — and 

n 

j/„(z) = I + (-l)"+ie*^+2^(-l)™e“^™, z>0. (7) 

m—1 

Here, it is understood that the nth pulse occurs at time 
tn = 5nt and 0 < < ... < (5„ < ... < ds < 1. In order 

to express the controlled decoherence function r„(t) in 
terms of its uncontrolled counterpart Toit), we simply 
relate to \yo{ojt)\'^ using iteration to find the 

following exact expression [^: 

n 

r„(t) = 2^(-i)"^+iro(t^) 

m—1 

n 

m—2 j<m 
n 

+ 2 ^ (-i)"^+"ro(t -+ (-i)"ro(t). (8) 

m—1 

In the case of one-sided noise, DD pulses are applied only 
to one qubit through the modified decay rate r(t). On 
the other hand, for two-sided noise, DD pulses are applied 
to both qubits. 

III. TIME-INVARIANT DISCORD 

Before starting to present our findings, let us first intro¬ 
duce some preliminary concepts and recall several known 
results form the literature, which have motivated us to 
investigate the time-invariant discord for dynamically de¬ 
coupled systems. The total amount of classical and quan¬ 
tum correlations in a bipartite quantum states can be 
quantified by the quantum mutual information, 

1{pab) = S{pa) + S{pb) — S{pab) 


(9) 
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where pab and pa(b) are the the density matrix of the 
total system and reduced density matrix of subsystem 
A(B) respectively and S{p) = —Tr{plog 2 p} is the von 
Neumann entropy. On the other hand, genuine quantum 
correlations of a more general type than entanglement 
are defined as quantum discord [40] 


'^{pab) = ^{pab) — C{pab), (10) 


where the classical correlations is given by |4T| 


CipAB) = max < S (pa) - V'pfc-5' {pA\k,) > (H) 

inn [ V J 

where {If®} is a positive operator valued measure 
(POVM) acting only on the subsystem B, and pA\k = 
TrB(Jl^PAB^k)/Pk is the remaining state of the sub¬ 
system A after obtaining the outcome k with probability 
Pk = TrABiA^PAB^k) i^ Ihe subsystem B. 

The curious phenomenon of frozen discord has been 
first realised in Ref. [5S] for a class of Bell-diagonal states 
which are locally interacting with Markovian dephasing 
environments. It was shown that quantum discord might 
remain constant at a non-zero value despite the effect of 
the noise for an initial period of time 0 < t t, during 
which only classical correlations decay. After the critical 
time t is reached, as discord is subject to noise-induced 
decay, classical correlations freeze. Consequently, this 
unanticipated phenomenon was named as sudden transi¬ 
tion between classical and quantum decoherence. On the 
other hand, for non-Markovian environments, the exis¬ 
tence of a transition time i crucially depends not only 
on the initial state of the bipartite system but also on 
the specific properties of the spectral density of the en¬ 
vironment |29| . Indeed, under certain conditions, it was 
demonstrated that a transition time t does not exist at all 
and therefore discord remains time-invariant throughout 
the whole time evolution [29]. 

As in other works studying time-invariant discord, we 
focus only on initial Bell-diagonal states of the form 


PAB = ^^|<I>±)($±| + ^i^|4/±)(vI-±| (12) 

where |4>) = (|00) ± \ll))/'/2 and \'i>) = (|01) ±\IQ))/V2 
are the four Bell states and |c| < 1. We consider two dif¬ 
ferent cases of local dephasing noise in presence of FDD 
in our work. Firstly, in case where both qubits A and 
B interact locally with identical dephasing environments 
having an ohmic type spectral density, the expressions for 
the mutual information T[pab) and classical correlations 


C{pab) are given by gH 
2 

C[pab] = E (13) 

f=i 
2 

nPABit)] = ^ l±i^log2[l + (-I)^C] 


i=i 


(14) 


where x(t) = max{e“^^(*\ c}. Secondly, we let only one 
of the qubits interact with the dephasing environment 
and assume that the other qubit is protected from de¬ 
coherence, in which case the expressions for the mutual 
information I[pab) and classical correlations C{pab) can 
be simply obtained by replacing with 

For the dephasing model we study, in case of un¬ 
controlled evolution, it is known that non-Markovianity 
manifests only when super-ohmic spectral densities are 
considered gSj. Note that the emergence of non- 
Markovian behavior here is directly connected with the 
decay rate 70 (t) being temporarily negative during the 
dynamics, which occurs for s > 2. At this point, we 
should also mention that although the approach we follow 
for characterizing non-Markovianity can be interpreted 
from several different angles based on the information 
exchange between the system and its environment, our 
results in fact do not depend on the specific choice of 
measure of non-Markovianity, since several well-known 
methods have been shown to qualitatively agree for the 
considered dephasing model mHH]. 

For simplicity, we consider the case of periodic DD 
pulses, applied at times = nAt with n = 1,2,3,.... 
We define two pulse intervals with respect to the onset 
of non-Markovian dynamics. In more detail, short pulse 
intervals are defined as At < f where i = tan(7r/s) and 
7 o(t) = 0 and conversely for long pulse intervals. Here, 
t denotes the first time instant after which the informa¬ 
tion flow is reversed for the uncontrolled dynamics. It 
was shown in Ref. |7] that optimal DD efficiency is al¬ 
ways achieved for Markovian regimes. In particular, for 
long pulse intervals, DD pulsing can lead to rapid deco¬ 
herence relative to the free dynamics. Indeed, if a pulse 
occurs during a period of re-coherence, initiated by the 
free non-Markovian dynamics, information flow will be 
inverted leading to decoherence. In the light of these re¬ 
sults, it becomes natural to relate the findings of Ref. 
[29] on the competition between non-Markovianity and 
DD to the creation of time-invariant discord. Specifi¬ 
cally, we intend to establish a set of necessary conditions 
for creating time-invariant discord. 

As continuously pulsed dynamics in the asymptotic 
time limit can not be defined analytically or numeri¬ 
cally, we focus rather on discord which remains “time- 
invariant” within a given experimental time rather than 
forever. Note that the values of invariant discord are 







4 



0 0.2 0,4 0.6 0.8 0 0.2 0.4 0.6 0.8 



0 0.2 0.4 0.6 0.8 0 0.2 0.4 0.6 0.8 


C c 

FIG. 1. (Color online) The black region shows the range of s 
and c parameters for which the discord is frozen forever for a 
free system evolving without pulsing. The gray region shows 
the range of s and c for dynamical decoupled systems with 
small interval spacing tUcAf = 0.3 (i) and long interval spacing 
tJcAf = 3 (ii). The plots in (a) are for the case where both 
qubits are affected by noise, as the plots in (b) are for single 
qubit noise. Outside these regions, one will always observe 
a transition from classical to quantum decoherence and thus 
no time-invariant discord. The final pulse is applied at = 
AmaxAt < 25a;J^ where Amax is the maximum number of 
pulses that can be applied within the time interval 0 < t < 

. Quantities plotted are dimensionless. 


given hy Q = {1 + c)log 2 (l -h c)/2 -h (1 - c)log 2 (l - c)/2, 
thus discord acquires a significant value for larger values 
of c. Fig. la (i) shows that, for two-sided local dephasing 
noise and short pulse intervals (e.g. WcAt = 0.3), time- 
invariant discord is created for a wider range of both c 
and s compared to the unperturbed case. Interestingly, 
we see that DD pulses allow for the existence of time- 
invariant discord also in the sub-Ohmic case (s < 1), 
for up to very high values of c. Generally, one can see 
from the plot that increasingly significant values of time- 
invariant discord (corresponding to large values of c) oc¬ 
cur for s < 2. Numerical investigation shows that these 
conclusions are independent of the specific pulse interval 
chosen, provided that At < t. 

In the asymptotic long time limit (tf ^ oo), based 
on numerical evidence, we conjecture that DD will only 
create invariant discord for s ;< 1. In more detail, as 
no degradation of the coherence is shown to occur within 
computable times, we assume that in the asymptotic long 
time limit (t/ —>■ oo), coherence is maintained for values 
close to unity. On the other hand, for s > 1, as t/ in¬ 
creases, the region of time-invariant discord will decrease 
as the coherence decays to smaller and smaller values. 

In Fig. la (ii), we show that time-invariant discord 
is always destroyed by DD techniques for large pulse in¬ 
tervals (e.g, WcAt = 3) in case of two-sided noise. Once 


again, we have checked numerically that this behaviour 
characterises the large pulse interval region. We conclude 
that for short pulse interval DD schemes and Markovian 
environments, specifically s < 1, one can create invariant 
discord for a larger range of initial states compared to 
relying on non-Markovianity alone as a resource. For the 
sake of comparison we show in Fig. lb (i) and Fig. lb 
(ii) exemplary plots illustrating the results of the same 
analysis when only one of the qubits is affected by de¬ 
phasing noise. Although we observe similar results to 
the previous case in general, we note that the regions of 
invariant discord, in both the uncontrolled and the con¬ 
trolled dynamics cases, become larger. However, we note 
that whereas the difference in the regions of s and c be¬ 
tween the one-sided and two-sided noise is significant in 
the case of unpulsed non-Markovian dynamics, the dif¬ 
ference is less pronounced when we implement the DD 
pulses. 

Let us now investigate the region of parameters c and s 
giving rise to time-invariant discord for intermediate val¬ 
ues of pulse interval WcAt. The behaviour of the bound¬ 
ary between the sudden transition and the time-invariant 
discord for different values of the pulse interval is dis¬ 
played in Fig. 2 (i) and Fig. 2 (ii) for the two-sided and 
the one-sided local dephasing noise, respectively. Anal¬ 
ogously to Fig. 1, we plot the transition boundary, be¬ 
low which, discord is time-invariant and classical corre¬ 
lations decay. Note that decoherence factors are set to 
unity initially at t = 0. Therefore, for a fixed s, if the 
initial state parameter 0 < c < 1 is chosen so that it 
is smaller than the minimum value of the decoherence 
factor at any point during the pulsed dynamics, 

there occurs no sudden transition and consequently dis¬ 
cord remains invariant in time at all times. However, 
when c > min(e“'"^*^) for a fixed Ohmicity parameter s, 
then this means that a sudden transition takes place and 
discord begins to decay after the critical time instant, 
which can obtained from the solution of the equation 
min(e“'"^*^) = c. Therefore, Fig. 2 actually demon¬ 
strates how the region of parameters s and c creating 
time-invariant discord are affected as the pulse interval 
changes between the short and long interval limits. Com¬ 
paring Fig. 1 and Fig. 2, we can see that when the pulse 
interval WcAt gets larger, the overlap between the regions 
of uncontrolled time-invariant discord and the controlled 
one becomes smaller. Finally, one can also observe that 
relying only on non-Markovianity as a resource of time- 
invariant discord becomes more preferable in general, as 
the pulse interval WcAt increases. 

We conclude our study by illustrating, for a specific ini¬ 
tial state, i.e. c = 0.5, how DD pulses affect the sudden 
transition and time-invariant discord dynamics for small 
and large pulse intervals with Markovian/non-Mar kovian 
regimes. Moreover, for completeness, we compare the 
behaviour of discord and classical correlations to the dy- 
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FIG. 2. (Color online) The boundary between the sudden transition and the time-invariant discord as a function of the Ohmicity 
parameter s and the minimum value of decoherence factor throughout the time evolution, for two-sided (i) and one-sided (ii) 
dephasing noise. We show pulse intervals tJcAf = 0.3 (red triangles), WcAt = 0.4 (blue squares), oJcAt = 0.5 (orange circles), 
tJcAf = 0.6 (green diamonds), oJcAf = 1 (blue plus signs), WcAt = 1.5 (purple crosses) and uJcAt = 2 (grey stars). The final 
pulse is applied at tf = NmaxAt < 25lo~^ where A^max is the maximum number of pulses that can be applied within the time 
interval 0 < t < 25a;ir^. Quantities plotted are dimensionless. 





FIG. 3. (Color online) Dynamics of entanglement (dotted 
green), classical correlations (solid red) and quantum discord 
(dashed blue) for i) short pulse spacing At = 0.3 and ii) 
long pulse spacing At = 3 (in units of uJc) and a) Markovian 
environments (s=1.01) and b) non-Markovian environments 
(s=4). The plots are for local dephasing noise on one qubit 
and c = 0.5. Quantities plotted are dimensionless. 

namics of entanglement as measured by concurrence: 

C = imax{0, |e“^^‘^(l —c)| —1 —c, |e“'"^*^(I-|-c)| — l-bc}, 

(15) 

where we assume that one of the qubits is locally in¬ 
teracting with the environment. For short pulse inter¬ 
vals (At < tan(7r/s)), and Markovian dynamics (s < 2), 


both classical correlations and entanglement are well pre¬ 
served while discord remains invariant at a non-zero value 
for the experimental time interval considered, as shown 
clearly in Fig. 3a (i). On the other hand, even if one of 
these conditions is not satisfied, time-invariant discord 
can vanish due to the sudden transition between classi¬ 
cal and quantum decoherence, and entanglement decays 
to zero. Quite intuitively, in order to preserve entan¬ 
glement through dynamical decoupling, a pulse needs to 
occur before entanglement decays to zero (see Fig. 3 (ii)). 


We note that, as one can see from Fig. 2, for c = 0.5 
time-invariant discord occurs only for smaller pulse spac- 
ings (wcAt ^ 0.6). Moreover, only for very small pulse 
spacing (wcAt « 0.3), can invariant discord be achieved 
for non-Markovian values (see Fig. 2 (i), red curve). 
Once discord become invariant, the value depends only 
on the initial state through the parameter c. 


Generally, in the Zeno regime and for divisible (Marko¬ 
vian) dynamics, dynamical decoupling protocols preserve 
quantum properties with increased efficiency for increas¬ 
ingly small pulse spacing [1S]-[1S]. The same holds 
for indivisible dynamics (non-Markovian) provided that 
t > At. In the short pulsed regime, it has been shown 
that as s increases, coherence is preserved less efficiently 
through DD protocols, revealing the detrimental influ¬ 
ence of non-Markovian dynamics [7j. For t < At, non- 
Markovian dynamics contribute to revivals in the dynam¬ 
ics and hence no general conclusion can be drawn. 




















6 


IV. CONCLUSIONS 


The persistence of quantum correlations in presence of 
noise induced by the environment is a remarkable phe¬ 
nomenon which may be crucial for those quantum infor¬ 
mation tasks using such correlations as a resource. In 
this paper we have studied the intriguing phenomenon 
of time-invariant discord with the aim of understanding 
whether the presence of DD pulses might lead to a wider 
range of parameters for which it occurs and if, in this 
way, one can obtain higher values of “protected” discord. 

We have performed an extensive analytical and numer¬ 
ical study exploiting an approach developed in Ref. [7] in 
a one-qubit scenario. While in the unpulsed regime time- 
invariant discord mainly occurs in the non-Markovian 
regime, our work proves that DD pulses allow for this 
phenomenon to take place also in the Markovian regime. 
In fact, our results show that, in the pulsed regime, the 
optimal conditions to achieve significant values of time- 
invariant discord persisting for the whole experimental 
time scales are obtained for Markovian reservoirs (s < 2) 
in the the short-pulse regime. Based on numerical anal¬ 
ysis, we conjecture that for discord to remain invariant 
to non-zero values forever, i.e. in the limit t/ —> oo, the 
Ohmcity parameter is further restricted to s < 1. On the 
contrary, large pulse spacing and free non-Markovian dy¬ 


namics are detriminental to the creation of time-invariant 
discord, in the pulsed scenario. 

Reservoir engineering is one of the most promising 
tools for fighting the loss of quantum properties, with 
prominent applications to realistic quantum technologies. 
Here we have considered two types of reservoir engineer¬ 
ing, one based on a modification of the reservoir spectral 
density (and hence the Ohmicity parameter), leading to a 
change in the Markovian character of the dynamics, and 
the second one exploiting pulses applied on the system 
(which, in turn, can also be seen as some sort of filtering 
of the spectral density). Our results show the optimal 
strategy to use, depending on the system-environment 
characteristic parameters. 

V. ACKNOWLEDGEMENTS 

C.A. acknowledges financial support from the EPSRC 
(UK) via the Doctoral Training Centre in Condensed 
Matter Physics under grant number EP/G03673X/1. 
G. K. is grateful to Sao Paulo Research Foundation 
(FAPESP) for the BEPE postdoctoral fellowship given 
under grant number 2014/20941-7. S. M. acknowledges 
financial support from the EU Collaborative project 
QuProCS (Grant Agreement 641277), the Magnus Ehrn- 
rooth Foundation, and the Academy of Finland (Project 
no. 287750). 


[1] H.-P. Breuer and F. Petruccione, The Theory of Open 
Quantum Systems, (Oxford Univ. Press, 2007). 

[2] L. Viola and S. Lloyd, Phys. Rev. A 58 , 2733 (1998). 

[3] L. Viola, E. Knill and S. Lloyd, Phys. Rev. Lett. 82 , 2417 
(1998). 

[41 K. Khodjasteh and D. A. Lidar, Phys. Rev. Lett, 95, 
180501 (2005). 

[5] K. Khodjasteh and D.A. Lidar, Phys. Rev. A 75 , 062310 
(2007). 

[6] T. E. Hodgson, L. Viola, and 1. DAmico, Phys. Rev. A 
81 , 062321 (2010). 

[7] C. Addis, F. Ciccarello, M. Cascio, G. M. Palma and S. 
Maniscalco, arXiv: 1502.02528 (2015). 

[8] A.Rivas, S. F. Huelga, and M. B. Plenio, Rep. Prog. 
Phys. 77 , 094001 (2014). 

[9] H. -P. Breuer, E. -M. Laine, J. Piilo, and B. Vacchini, 
arXiv:1505.01385 (2015). 

[10] B. Bellomo, R. Lo Franco, G. Gompagno, Phys. Rev. 
Lett. 99, 160502 (2007). 

[11] A. Rivas, S. F. Huelga, M. B. Plenio, Phys. Rev. Lett. 
105 , 050403 (2010). 

[12] H.-P. Breuer, E.-M. Laine, J. Piilo, Phys. Rev. Lett. 103 , 
210401 (2009). 

[13] S. Luo, S. Fu, and H. Song, Phys. Rev. A 86, 044101 

( 2012 ). 

[14] B. Bylicka, D. Ghruscihski, S. Maniscalco, Sci. Rep. 4 , 
5720 (2014). 

[15] S. C. Hon, X. X. Yi, S. X. Yu, and G. H. Oh, Phys. Rev. 
A 83, 062115 (2011). 


]16] X.-M. Lu, X. Wang, and C. P. Sun, Phys. Rev. A 82, 
042103 (2010). 

[17] F. F. Fanchini, G. Karpat, B. Qakmak, L. K. Gastelano, 
G. H. Aguilar, O. J. Farias, S. P. Walborn, P. H. Souto 
Riberio, and M. G. de Oliveira. Phys. Rev. Lett. 112 , 
210402 (2014). 

[18] S. Haseli, G. Karpat, S. Salimi, A.S. Khorashad, F. F. 
Fanchini, B. Qakmak, G. H. Aguilar, S. P. Walborn, and 
P. H. Souto Ribeiro, Phys. Rev. A. 90 , 052118 (2014). 

[19] D. Ghruscihski and S. Maniscalco, Phys. Rev. Lett. 112 , 
120404 (2014). 

[20] F. F. Fanchini, G. Karpat, L. K. Gastelano, and D. Z. 
Rossatto, Phys. Rev. A 88, 012105 (2013). 

[21] C. Addis, B. Bylicka, D. Ghruscihski, and S. Maniscalco, 
Phys. Rev. A 90 , 052103 (2014). 

[22] R. Horodecki, P. Horodecki, M. Horodecki, and K. 
Horodecki, Rev. Mod. Phys. 81 , 865 (2009). 

[23] K. Modi, A. Brodutch, H. Gable, T. Paterek, and V. 
Vedral, Rev. Mod. Phys. 84, 1655 (2012). 

[24] V. Madhok and A. Datta. Phys. Rev. A 83 ,032323 
(2011); D. Gavalcanti, L. Aolita, S. Boixo, K. Modi, M. 
Piani and A. Winter, Phys. Rev. A. 83 , 032324 (2011); 
A. Streltsov, H. Kampermann and D. Bruss, Phys. Rev. 
Lett. 108 , 250501 (2012); T. K. Chuan, ,H. Maillard, K. 
Modo, T. Paterek, M. Paternostro and M. Pianai, Phys. 
Rev. Lett. 109 , 070501 (2012). 

[25] S. Boixo, L. Aolita, D. Cavalcanti, K. Modi, A. Winter, 
Int. J Quant. Inf. 9, 1643 (2011). 



7 


[26] M. F. Cornelio, M. C. de Oliveira and F. F. Fanchni, 
Phys. Rev. Lett. 107, 020502 (2011). 

[27] A. Streltsov, H. Kampermann and D. Bruss, Phys. Rev. 
Lett. 106, 160401 (2011); T. Werlang, C. Trippe, G. A. 
P. Ribeiro and G. Rigolin, Phys. Rev. Lett. 105, 095702 
( 2010 ); 

[28] L. Mazzola, J. Piilo and S. Maniscalco, Phys. Rev. Lett 
104 200401 (2010). 

[29] P. Haikka, T. H. Johnson, and S. Maniscalco, Phys. Rev. 
A 87, 010103(R) (2013). 

[30] L. Mazzola, J. Piilo, and S. Maniscalco, Int. J. Quant. 
Inf., 09 981-991 (2011). 

[31] R. Lo Franco, A. D’Arrigo, G. Falci, G. Compagno, E. 
Paladino, Phys. Rev. B 90, 054304 (2014). 

[32] S. Singha Roy, T. S. Mahesh, and G. S. Agarwal, Phys. 
Rev. A 83, 062326 (2011). 

[33] F. F. Fanchini, E. F. de Lima, L. K. Castelano, Phys. 
Rev. A 86, 052310 (2012). 

[34] A Rosario et al 2012 J. Phys. B: At. Mol. Opt. Phys. 45 
095501. 


[35] Jian-Song Zhang and Ai-Xi GHen, J. Phys. B: At. Mol. 
Opt. Phys. 47 21 (2014). 

[36] J. Luczka, Physica A 167, 919 (1990). 

[37] G. M. Palma, K.-A. Suominen and A. K. Ekert, Proc. 
Roy. Soc. Lond. A 452 567 (1996). 

[38] J. H. Reinga, L. Quiroga and N. F. Johnson, Phys. Rev. 
A. 65, 032326 (2002). 

[39] G. S. Uhrig, Phys. Rev. Lett. 98, 100504 (2007); G. S. 
Uhrig, New J. Phys. 10, 083024 (2008). 215505. 

[40] H. Ollivier and W. H. Zurek, Phys. Rev. Lett. 88, 017901 

( 2001 ). 

[41] L. Henderson, and V. Vedral, J. Phys. A 34, 6899 (2001). 

[42] S. Luo. Phys. Rev. A 77, 042303 (2008). 

[43] P. Facchi, S. Tasaki, S. Pascazio, H. Nakazato, A. Tokuse 
and D. A. Lidar, Phys. Rev. A 71 022302 (2005). 

[44] P. Facchi, H. Nakazato and S. Pascazio, Phys. Rev. Lett. 
86 2699 (2001). 

[45] P. Facchi, D. A. Lidar and S. Pascazio, Phys. Rev. A 69 
032315 (2004). 



